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I. Introduction 



The production of charmonium pairs in the same reaction is expected to be exceedingly rare in QCD. However, the 
NA3 collaboration has measured ip^p pair production in multi-muon events in 7r~Pt interactions at 150 and 280 GeV 
and in 400 GeV pPt interactions The integrated ip'ip production cross section per nucleon, cr^.^, is 18 ±8 pb at 
150 GeV and 30 ± 10 pb at 280 GeV for pion-induced and 27 ± 10 pb proton-induced interactions at 400 GeV. This 
cross section is reduced relative to the single J/tp production cross section, cr^, by a^^/a^ « (3 ± 1) x 10"'' in the 
pion-induced events. 

The fraction of the projectile momentum, = pf^/pheiam carried by the ijiij) pair is always very large, > 0.6 
at 150 GeV and x^^ > 0.4 at 280 GeV. Additionally, single J /if) mesons in the VV' pairs have > 0.15. When 
the data from both pion beams are combined, (cc^^) = 0.66 and {x^) = 0.33. As we shall show, perturbative gg 
and qq fusion processes produce central ij^ip pairs with lower average momentum fractions^. The average invariant 
mass of the pair, (M^^) = 7.4 GeV, is well above the -0-0 production threshold, 2m^ = 6.2 GeV. In fact, all of the 
observed events have an invariant mass greater than 6.7 GeV. Additionally, the average transverse momentum of the 
pair is quite small, PT.^^p = 0.9 ± 0.1 GeV, suggesting that pair production is highly correlated, occurring in the 
same interaction. The pN events have a somewhat lower invariant mass, (M^^) w 6.8 GeV, with several events near 
threshold. The x^p^p distribution was not reported. 

In this paper, we shall investigate a number of leading-twist and higher-twist ipij} pair production models. Several 
leading- twist production mechanisms have been discussed previously, including 0(q!^) i/'V' production via gg fusion 
and qq annihilation (|-||], BB production and decay, BB ip^X and the production of a 2'^~^cccc resonance 
that decays into ipip pairs 0. We review double J/ip hadroproduction by these leading-twist mechanisms, updating 
the results with more recent parameters and nonscaling parton distributions. These models are unable to produce 
ip'ip pairs with large x^^. 

We then turn to a higher-twist mechanism that easily produces fast QQ pairs: intrinsic heavy quark components 
in the projectile wavefunction In a companion paper, Ref. we evaluated the probability for two J/tp^s to 
originate from an intrinsic \udcccc) Fock state. This model was successful in reproducing the general features of the 
data. In this paper, we also consider several other sources of double J/ip production by intrinsic heavy quark states. 
An intrinsic bb pair, \udbb)^ could decay into "00 pairs, as in leading-twist bb production. Additionally, an intrinsic 
\udccbb) state could produce tpip pairs from a -B — > J/tpX decay combined with the coalescence of the cc pair into a 
J/ip. Finally, we calculate ^p^p pair production through leading-twist cc production with the projectile in an intrinsic 
cc state. Examples of all the processes considered are shown in Fig. 1. 

In addition to estimating the ipip production cross section for all of the leading-twist and higher-twist mechanisms, 
we compare our results to the combined NA3 pion-induced data and the invariant mass distribution from the 400 
GeV proton-induced production. We also discuss ipip pair production by 800 GeV proton beams, measurable in 
ongoing fixed-target experiments. To facilitate comparison with the data, we present all the momentum distributions 
in the laboratory frame and normalize the pion-induced calculations to the combined 150 and 280 GeV data and the 
proton-induced calculations to the 400 GeV data. The -00 pair mass and momentum distributions are normalized to 
the pair rate while the J/0 momentum distributions are normalized to the single rate, twice the pair rate. 



II. Leading-Tv^fist Production Mechanisms 

In this section we shall review the leading-twist QCD predictions for double quarkonium hadroproduction, O(a^) 
ip'ip production and bb production and decay. We also estimate double J/tp production assuming the existence of a 
2+"'"cccc resonance. These processes are sensitive to the parton distributions in the pion and nucleon. We use two 
sets of recent leading-order parton distribution functions in our calculations, GRV LO, available for both nucleons 
pof and pions |]ll|, and Duke-Owens 1.1 |l^ for the nucleon with Owens I for the pion. The results we present 
are calculated for a proton target, but at these energies and with these scale choices, the 7T~p cross sections can be 
« 50% larger than the corresponding N cross sections since the contribution from pion valence quark annihilation 
is important. There is no significant difference between the pp and pN production cross sections. 

Since the calculations are leading order, a K factor is usually invoked to account for next-to-leading order corrections. 
For example, in bb production, K = a{al)/a{al) « 1.2 — 2.5, depending on the choice of mi, and the scale parameter 



'^Single J/0's produced at rest in the tt~ N center of mass have laboratory momentum fractions of = 0.18 at 150 GeV and 
0.14 at 280 GeV. Likewise, assuming M^^ = 7 GeV leads to laboratory fractions x^^ = 0.42 at 150 GeV and 0.3 at 280 GeV 
for 0i/) pairs at rest in the n'N center of mass. 
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fib 0. Indeed, for fib = rub, a{mb = 4.5 GeV)/cr(m6 = 5 GeV) « 3.7 at 0{al) while for rUb = 4.75 GeV, a{fib = 
0.5nib)/<7{fib = ^nib) « 6 at 0{al) and 3 at 0{al). Thus over a reasonable range of rUb and fib, the leading order bb 
production cross section can vary by an order of magnitude at the energies studied here. 

There is thus considerable uncertainty in the normalization of the leading-twist production cross section. In 
our analysis we will quote a representative value of a^^,. The most important test of these models, however, is the 
shape of the momentum and invariant mass distributions since these are not strongly affected by the target type or 
the scale parameters. We will find that while some of these mechanisms can account for the magnitude of cr^^ within 
a factor of 2-3, none of them can fully describe the characteristics of the momentum and invariant mass distributions. 



1. Production at C'(af) 

Leading- twist double J/^p production has been calculated to leading order in as, 0(0;^), in Ref. (^1). In these 
caclculations, the J/i/j's are treated as nonrelativistic bound states with the appropriate spin and angular momentum 
quantum numbers, as in the color-singlet model of quarkonium production [ p^ . Examples of gg — *■ Tpip production 
diagrams are shown in Fig. 1(a). There are 36 diagrams in all, some vanish due to color conservation and others can 
be related by crossing [El. The qq — s- ^pip diagrams |^,^ are shown in Fig. 1(b). After four- momentum conservation is 
accounted for, the 0{aY) production cross section may be written as 

^-^^^ = E(^.(-)^.(-) + ^.(^i)^,(^.))f I,. + ^.(xOF.(^.)f I.. , (1) 

where cci = mT,^(e^i -t-e''^)/y's, X2 = mT,^{e~^^ +e~y'^) / ^ andm^^ = m^+p^. The parton momentum distribution 
functions, F{x) — xG{x), are evolved to — 4m^ so that when the renormalization and factorization scales are 
assumed to be equal, « 0.21. 

The subprocess cross sections for xj:^ production to order 0{a^) from qq annihilation and gg fusion are Q| 

da, 2 /167ra^\'^ |4'(0)n*(0)P 1 9 9 4, 

dt 37r V 3 / 

^1 _ 3 /^g.y |vl/(0)nvl/(0)P s^TWCniz) 

dV'' 47r V 3 / (£ - )4(w - )4 ' ^> 

where y — {tA~ u)/2s, z = [t — u)/2s and s = XaXbS, t — —ximT,ipy/se~y^ + m^, and u = — X2mT. 4,^/3 ey'^ + m^. The 
J/ip wavefmiction at the origin, |5'(0)p, is fixed by the electromagnetic decay width, V{J/ip e+e~). We use the 
most recent determination, |*(0)p = 0.0689 GeV^ Q. The polynomials C„(z), from Ref. [Q, are 

Co = I (39813122^2 ^ 2580482^" - 1186048z^ + SbbSAz^ 

+ 98U4z'^ - 13568^2 + 335) 
Ci = ^(1935360z^° - 596224.2^ - 403840z^ + 168084z^ 

- 14288^2 + 137) 

C2 = -(1134336^^ - 1010944z^ -I- 239008z^ - 11472^^ - 1) 
3 

64 

C3 = — (-10176z^ + 688z'' + 828^2 -f- 5) 
C4 26(16zV8z2 + 1) . 

The predicted cross sections for ir^p — > ipip are 7.4 pb at 150 GeV and 16.5 pb at 280 GeV, « 30-50% of the 
measured total. At 400 GeV, the predicted pp production cross section is 7.2 pb, « 30% below the data, 

increasing to 26 pb at 800 GeV 0. 



^In Ref. 1^, a smaller value of the wavefunction, |\l/(0)|^ = 0.0387 GeV^, and a larger value of the strong coupling constant, 
Oa — 0.32, were used, leading to similar conclusions about the magnitude of the production cross section. 
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The calculations are compared with the 7r~iV — > ipiJjX data in Fig. 2. The results are not strongly dependent on the 
choice of parton distributions. The momentum distributions are relatively broad, with (x^) = 0.26 and (x^^) = 0.5 
at 150 GeV; (x^) — 0.21 and {x^^) = 0.4 at 280 GeV. Note that {x^^) is significantly less than required by the ^ip 
data. In fact, all the measured pairs have a larger x^^ than these averages. Also 20% of the single J/ip^s in the pairs 
have x^ > 0.5, larger than predicted by leading order QCD ipip production. The pp -tpip calculations, shown in Fig. 
3, predict narrower distributions than the ir^p calculations: (x^) — 0.15 and (x^^,) = 0.3 at 400 GeV; {x^) — 0.11 
and {x^ip) = 0.24 at 800 GeV. We find (M^^) — 7 GeV for both projectiles. This is in good agreement with the 400 
GeV proton- induced data where the pair mass is low. However, the ir^ N data suggests that production occurs above 
threshold. Thus although the absolute cross sections are not significantly below the data, the mass and longitudinal 
momentum distributions from this model fail to describe the data. 



2. Leading-Twist hb Production 

We estimate the ip^f) pair production cross section from leading- twist bb production and decay. At ©(a^), bb 
production proceeds by gg fusion and qq annihilation, as shown in Fig. 1 (c) . The bb production cross section has the 
same form as Eq. (1) but the subprocess cross sections are ||l7[| 

da ^ Trgg cosh(yb - j/^) + mg/mg 
di^'" 9mt (l + cosh(yb-y^))' ' 



da _ iral 8 cosh(yb - y^) - 1 (_ ^ , ^ , 2m^ 2m'' 

dl'^" " 96m4 (1 + cosh(yb - y^)f 



cosh(yb - J/.) + ^ - ^ . (5) 



The parton distributions are evaluated at — 4to^ with to& = 4.75 GeV. The predicted rates of open bb production 
are a^^ij^-p) = 230 pb at 150 GeV and 2 nb at 280 GeV. We also find a^^{pp) = 250 pb at 400 GeV and 1.3 nb at 

800 GeV. As previously discussed, changing m.\, and Q ~ \ih can change the predicted (t*2* by a factor of five to ten. 

Flavor excitation has been considered as a source of b quark production in Ref . |^ . These quark-gluon scattering 
graphs are considered part of the K factor which increases the cross section K k.2 |i^, as already discussed. In 
fact, if the uncertainties are varied to produce the maximum production rate, the 'n~ p production cross section could 
be increased to ~ 10 nb, not the 100 nb predicted in Ref. A later study ruled out 66 production as a significant 
source of '0V' P^ir production due to the small production cross section, further reduced by the BB — + '\\)'\\} decay. 
Indeed, including the square of the inclusive branching ratio, BR{B J/ijjX) w 1.3%, |l^ the results are a factor 
of 10 to 100 below the data^. Our predicted 7r~p ipi' production cross sections are then 0.028 pb at 150 GeV and 
0.24 pb at 280 GeV while the corresponding pp ipip cross sections are 0.03 pb at 400 GeV and 0.16 pb at 800 GeV. 

Measurements of charm hadroproduction show that the longitudinal momentum distribution of final-state D mesons 
are similar to the predicted charmed quark distributions. Using a fragmentation function that reproduces D production 
in e~^e~ annihilation JlOt softens the D distribution further. However, in hadroproduction, comoving light partons can 
combine with the charmed quarks without significant momentum loss by the charmed quark so that the fragmentation 
can be described by a delta function, D£,/c{z) w 5{z — 1), in matter pOf . Since the b quark fragmentation function 
from e"'"e~ annihilation predicts less momentum loss |2l| ], the delta function approximation of fragmentation should 
be even better for b quarks. Therefore x^g — x^j^ and xb — Xb- However, the mass distributions are kinematically 
affected since the B meson mass reduces the available phase space for high mass BB pairs. We find that the average 
BB mass is approximately 1 GeV larger than the average bb mass and that both increase with incident energy. This 
predicts a larger ^pip invariant mass than 0{q.'^) ipip production since each J is produced by an isotropic decay, 
leading to a larger rapidity gap between the J/i/j's and thus higher i/'V' pa-ir masses. Note that the ■0V' pair masses are 
larger if the b quarks are assumed to decay into J/ip^s directly. However, the single and pair J/iJj distributions are 
only weakly dependent on whether the parent particle was a b quark or a B meson. The x^^ and x^^ J ftp distributions 
have larger average values and narrower widths than the 0{a'^) calculations. 



^We assume that all b quarks produce final states that can decay to J/ip's to get an upper bound on the production cross 
section and study inclusive B decays. Choosing an exclusive decay channel, such as B ^ J/tpK, reduces the branching ratio. 
Other assumptions of how the final-state hadrons are produced would also reduce the J/tp rate. 
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We calculate the B J/ipX decay distributions by a Monte Carlo using J/tp momentum distributions calibrated 
to the B decays measured at CLEO The J/'0 and ip^p distributions are not smooth due to the finite statistics of 
the Monte Carlo. 

The B and J/tp laboratory frame distributions are given in Figs. 4 and 5. The B and BB longitudinal momentum 
distributions are calculated with both sets of parton distribution functions. The invariant mass, J/ip and tpip momen- 
tum distributions are shown for the GRV LO partons only to illustrate the differences arising from assuming either 
that the b quarks decay directly to J/iJj or that the quarks first hadronize and then decay. The average momentum 
fractions and invariant masses of the B mesons and J/V''s calculated with the GRV LO parton distributions are shown 
in Table I. In addition to the shape differences between the calculations and the data, the predicted normalization is 
a factor of « 100 smaller than required by the data. Thus we can conclude that leading-twist bb production is not an 
important source of ipip production at these energies. 



3. 2^~^cccc Resonance Production 



In the resonance model of Ref . , ip^p pairs with M.^^ = 7 GeV are produced by a Drell-Yan type mechanism where 
two gluons fuse in a color vector-meson dominance model. Each gluon couples to a cc pair in the color-octet-vector 
representation. The cc pairs then couple to a 2"'""'"cccc state which subsequently decays to V'V' pairs. The basic process, 
which has been generalized to a number of vector- meson final states, is shown in Fig. 1(d). We follow Ref. H in our 
calculations. The ipijj pair longitudinal momentum distribution in the laboratory frame is 



da 



1 

SPL 



(6) 



where Gg{x) represent the gluon number densities, z is the fraction of the laboratory momentum, p^, carried by the 
'0V' pair, W is the center of mass energy of the two gluons, and W^,jjj = 4m^. The gluon fusion cross section is 
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64 US-kW 3 



7 / 2fc2 

1 + 7^ + 



3m2 



2fc4 
15to^ 



7^ h"^ 



(7) 



where a^^ = a/\/3 is the resonance decay constant, b^^ = y^2/5asa{ATr/ f^)'^ is the coupling between the gluons and 

the 2++ state, a ~ -\/30, and (47r//^)^ ~ 0.02, the color vector dominance constant. The width of the -00 resonance 
is 




(8) 



where k = O.h^^W^ — 4m^ is the J/t/j momentum in the center of mass. A single J/'ip in the ipij} pair has the 
laboratory momentum distribution 



da 

dp^ 



1 



dW 



dz- 



2yr 



^z'^ + W'^Ip\ 



(9) 



where vp = zj ^ z"^ + W^/pj^ and the limits on z are determined from the angle, 6, between the J/i^^s, —1 < 
W{p^ — zpL/2)/{ky/z^pj~rW^) < 1. The angular distribution resulting from the decay of the 2++ state is 



dar, 



1 



7 



fc2 



dcost 
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+ ■ 



5 4 2/,, 4 

- - - COS 9 + cos ( 



5 cos^ 



(10) 



7^ /i2 



{W-M^^f+Tl^{W)/4 



In the original model of Ref. ||^, scale invariant parton distribution functions with as = 0.18 were used. Repeating 
the calculation using the same gluon distributions, we find agreement with their results, within one standard deviation 
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of the measured cross section^. However, if the GRV LO gluon distributions are used at Q = M^^ = 7 GeV, the 
cross sections are reduced significantly. We find 1.7 pb at 150 GeV and 10 pb at 280 GeV for pion-induced production; 
12.8 pb at 400 GeV and 55.2 pb at 800 GeV from proton-induced production. The momentum and mass distributions 
are shown in Figs. 6 and 7. Although the model predicts narrower momentum distributions than the 0{a^) process, 
the averages are similar. However, neither mechanism produces fast ^if) pairs. 



III. Higher-Twist Production: Intrinsic Heavy Quark Mechanisms 

In the following sections, we will discuss the production of one or two J/ip^s from intrinsic heavy quark states in the 
projectile wavefunction. The wavefunction of a hadron in QCD can be represented as a superposition of Fock states 
of quarks and gluons, e.g. the 7r~ wavefunction includes \ud), \udg), \udQQ), etc. components. (Since our calculations 
are applied to pion and proton projectiles, we use the notation \nvQQ) for the heavy quark Fock state where ny = ud 
for a TT" and uud for a proton.) When the projectile scatters in the target, the coherence of the Fock components is 
broken and its fiuctuations can hadronize |p,p3|. For example, intrinsic cc fluctuations can be liberated, provided the 
system is probed during the characteristic time, At — 2piab/M^-, that such fluctuations exist. 

Microscopically, the intrinsic heavy quark Fock components are generated by virtual interactions such as gg — > QQ 
where the gluons couple to two or more projectile valence quarks. The probability to produce QQ fluctuations scales 
as a^g{rriQ) / mQ relative to leading- twist production [2^ . Thus intrinsic heavy quark production is higher twist. 

The dominant Fock configurations are not far off shell and thus have minimal invariant mass, JVP — rn^^/xi, 
where i = "^f transverse mass of the i^^ particle in the configuration. Intrinsic QQ Fock components 

correspond to configurations with equal rapidity constituents. Thus, unlike sea quarks generated from a single parton, 
intrinsic heavy quarks carry a large fraction of the parent momentum ||]. In fact, if the intrinsic QQ coalesces into a 
quarkonium state, the momentum of the two heavy quarks is combined and this final state carries a dominant fraction 
of the projectile momentum. 

There is substantial circumstantial evidence for intrinsic cc states. For example, the charm structure function of 
the proton measured by EMC is significantly larger than predicted by photon-gluon fusion at large x bj [P5| . Leading 
charm production in nN and hyperon-iV collisions also requires a charm source beyond leading twist ]26| , |27[ | . The 
NA3 collaboration has also shown that single J/ tp production at large xp is greater than expected from leading- twist 
production [2^ . 

The probability distribution for a general n-particle heavy quark Fock state as a function of x and kr is written as 

Ylt.dx.d-'kT, ''"''^^'''qq> K-e:=i("^I,.M))^ ' ^ ' 

where Nn, assumed to be slowly varying, normalizes the probability. The assumed constituent quark masses are 
niq — 0.3 GeV, mc = 1.5 GeV, and mb = 4.75 GeV. While the /ct dependence is needed to calculate the mass 
distributions, an average k\ can be used to calculate the x distributions. Thus, eq. (11) can be simplified to 

where rhf = m|-f (fc^ ,j) is the effective transverse mass. Assuming (fc^ ,j) is proportional to the square of the constituent 
quark mass, we adopt the effective values fhc = 1.8 GeV and fhq — 0.45 GeV, as in our previous work |^,|9) and 
use rhb = 5 GeV. Note that the resulting x distributions are independent of the Lorentz frame. We will show results 
involving B decays in the laboratory frame. Eqs. (11) and (12) can be generalized to an arbitrary number of light 
and heavy partons in the Fock state so that the probability distribution for the production of e.g. two heavy quark 
pairs is, 

^/^Q-Q- = N^atiM^ ^ )at{M^ ^ ■ (13) 

YYl^^dx^d^kT., -Er=lHT^M))' 



*The gluon convolution was assumed to be G'll{x\)G'^'^{x2) + Gg^ (a;i)Ggj^ (2:2), twice as large as usual For consistency 
with Ref. iM], we included this factor in our calculations. However, this difference should be noted. 
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The intrinsic QQ production cross section can be related to the probabihty, PiQ, by 

aiQ(/i7V)=FiQaj;V^ . (14) 

In their single J measurements, the NA3 collaboration separated the nuclear dependence into a "hard" contribution 
with a nearly linear A dependence at low xf and a contribution scaling as for nA interactions, characteristic 

of soft interactions, called "diffractive" . The soft interaction scale parameter, /i^ ~ 0.2 GeV'^, was fixed by the 
assumption that the diffractive fraction of the total production cross section measured by NA3 |2^ is the same for 
charmonium and charmed hadrons. Therefore, we obtained (Tic(7rA^) « 0.5 /ib at 200 GeV and aic{pN) « 0.7 /ib |2^] 
taking Pjc = 0.3% from an analysis of the EMC charm structure function measurements p^. 



1. Production Prom \nvbb) Pock States 

We first examine tpi/j production from a \nvbb) configuration. Because of the scaling properties of the intrinsic QQ 
production probability |27| ] the intrinsic bb production cross section is 

a,AhN)=aUhN)(^!^f^y . (15) 



\rhb as{M, 
Thus aib{hN) « 2-3 nb and 

af^^'f' {tt- N) = [BR{B ~> J/i/^X)] Vib(7r" A^) « 0.36 pb . (16) 

The corresponding proton-induced cross section, cr^{pN) w 0.51 pb, is larger. These results are larger than our 
quoted leading-twist bb production cross sections. However, this difference is strongly dependent on the leading-twist 
parameters. Also, the intrinsic bb cross section is proportional to the hN inelastic cross section and changes slowly 
with ^/s while the leading-twist cross section has a steep ^/s dependence, particularly near threshold. Although bb 
production and decay is frame independent, we compare the resulting J/^/^'s with the data in the laboratory frame. 

There are two ways of producing B mesons from the intrinsic bb pairs. The first is by standard fragmentation 
processes. We will assume that the momentum of the quark lost through fragmentation is small so that a delta 
function can be used as the fragmentation function, as discussed for leading-twist production. This assumption works 
well for nonleading D meson production |^,^. Then the meson and quark distributions are identical. The quark 
can also coalesce with a projectile valence spectator to produce mesons if the projectile has the corresponding valence 
quark. The coalescence mechanism introduces flavor correlations between the projectile and the final-state hadrons 
and produces B^s with a larger fraction of the projectile momentum^. Such a coalescence model has been used to 
successfully describe production asymmetries We denote B mesons produced by fragmentation with 

Bp and those produced by coalescence with Be- The production processes are illustrated in Fig. 1(e). 

If we assume that the b quark fragments into a B meson, the B distribution is 

diF^' f dP' 

, ' = I dzY\dxi- — DB/b{z)5{xBF - zxb) , (17) 

dXBp J fJi dxi . . . dXn 

where n = 4,5 for pion and proton projectiles and Dg/b^z) = S{z — 1) is the b quark fragmentation function. The 
fragmentation mechanism produces B mesons with 25-30% of the projectile momentum. While this is comparable to 
or larger than the laboratory fractions carried by B's produced by leading- twist fusion, it is not as large as the result 
for Be production. The Be distribution is 

dP' f df^' 

, ' = Y\dxi- — 6{xbc - Xb - xi) . (18) 

dXBc J , dxi...dxn 



^In this model, with a tt projectile, B (ub) and B^{db) can be produced by both coalescence and fragmentation while 
B^{ub)/B (dh) can only be produced by fragmentation. In a proton projectile, B^/B'^ are produced by coalescence and 
fragmentation. We have not made any assumptions of the relative B meson production rate by coalescence and fragmentation. 
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In this case, the Be can carry 40-50 % of the projectile momentum. The increase over Bp production is due to the 
Ught valence quark. The B and J/V' distributions are shown in Fig. 8(c) and 8(d) for the pion and Fig. 9(c) and 9(d) 
for the proton. 

The ^l^if) pairs can be produced either from the decays of the BE states produced by fragmentation of the intrinsic 
bb pair, the combination of fragmentation and coalescence, or, from a pion projectile only, double coalescence. The 
BpBp distribution is 

* J dzbdz-^JJdxi DB/b{zb)D^^-^{z-^)S{xBp - ZbXb)S{xg^ - z-^x^) . (19) 



These tpip pairs have the smallest average momentum fractions of all tpip pairs produced from the intrinsic \nvbb) 
configuration. If the final-state pair is BcBp, 



/ dP' 
dz J\ dxi D^^^{z)5{xbc -Xb- xi)5{x^^ - zx^) . (20) 



An initial \uudbb) configuration provides no valence antiquark for Be production, thus the tpi/j pairs resulting from 
Be Bp decays carry the largest fraction of the proton momentum. However, when the projectile is a pion, such double 
coalescence is possible. All the pion momentum is transferred to the BcBc pair 



j n dxt^'^dxi ^^^'^" ^^''^^^Bc ~^b- X2) , (21) 



^^BcBc ■> i=l 

i.e. (a;g^,-g^,) = 1. The Be Be ipipX decay results in high momentum ^tp pairs. The BB distributions are shown 
in Figs. 8(a) and 9(a). 

The BB mass distributions can determine the relative probability for the hadronization mechanism since the 
available phase space depends on the number of light quarks involved. The mass distribution of the bb quark pair is 



^ = / n dx^d'kT, ^rf^fc^,,, ^.^^^^^^^,^^ / (x,, -Xb- X,) (22) 



X S{kT,b + kj,-^-kj,f;^)S 



'M, 



2 



'T,bb '^T,b 



Xb 




We find {M^-^ = 11.75 GeV. The corresponding M^^ distributions are broad with (M^^) = 8.0 GeV. A significant 
part of the hb mass distribution lies below the BB production threshold, 10.56 GeV. In order to compute BpBp 
production through delta function fragmentation, we restrict mx^b to be larger than vriB- The mass cut reduces the 
yield by a factor of ~ 20 and increases the average mass to {^Bf'Bp) ~ ■'^^■^ ^^''^ both projectiles. This is the 
largest average BB mass because although a mass threshold has been imposed, the light quark transverse momentum 
has not been affected. The avciragc; i/'i/' mass is reduced to 7.65 GeV, independent of the projectile. 
When Be Bp pairs are created, one of the projectile quarks is included in the heavy meson pair, 

d^k^ BcsA'^BF - 5.28GeV) (23) 

5{kT,l + kxfi - kT,Bc)^{XBc - Xb- Xi) 

J^b -i^l I 

T,Bf ' ^ 




^BcBf 



S{kT,Bc + kj, B^ - k^^BcBF^^^^BcBF ~ ^^^c - XBf) , 



reducing the invariant mass of the BB to 12.2 GeV for pion and 12.4 GeV for proton projectiles. When the projectile 
is a pion, the Be Bp pair momentum has to be balanced by the remaining valence quark, resulting in a stronger 
correlation than in a proton with two left over valence quarks and a smaller pair invariant mass. The V'V' Pair mass 
is also reduced to 7.5 GeV for a pion projectile and 7.6 GeV for the proton. 
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The strongest momentum correlations arise from BcBq production since all the pion momentum is transferred to 
the meson pair, i.e. x^^-g^ = 1 and kj, ^^-^^ = 0. The mass distribution for these pairs can be simplified to 

— — = / dxidx2 '^^^ -dm%^drrT^ (fkT^bd^kr^-^(fkj,-g^, (24) 

"-'"o -p J XBc\'- XBc) 



X 5 



'^T.Bc 






XBc 




Xl j 




2 


9 

m — 

T,Bc 


XBc 


^-^Bc 



1-XBc '^-XBc~X2 X2 




.2 



Xl X2 XBc—Xi I — XBC—X2J 

The X integrations are done in the delta functions, leaving the fc^ and B meson mass integrals to balance the pion 
invariant mass, producing the steepest pair mass distribution with (Mg^-^^) — 11.8 GeV. The tpTp mass is also 
reduced to (M^^) = 7.05 GeV. The strong phase space restrictions decrease the relative production probability for 
this mechanism, making double coalescence rare. However, the detection of such states would be a dramatic proof of 
the existence of the coalescence mechanism. 

With these pair mass distributions, the average tptp momentum is decreased by 5-10% with respect to the average 
obtained assuming that the b quarks decay to J/V'- The resulting average momentum fractions and invariant masses 
are given in Table II. 

Although the intrinsic \nvbb) configuration produces fast ^''4' pairs, with 45-80% of the pion momentum, the ^p'lp 
production cross section is nearly a factor of 100 too small to account for the NA3 data, primarily due to the small 
branching ratio for B ^ J / tpX decay. 



2. Production From Doubly Intrinsic QQ Configurations 

It is possible for configurations with more than one intrinsic heavy quark pair to exist in the projectile, shown in 
Fig. 1(f). Even though the probability for such configurations is reduced relative to those with only one intrinsic QQ, 
the ipip production cross section could be larger from \nvcccc) or \nvccbb) configurations than the \nybb) because the 
square of the B decay branching ratio does not enter. We now discuss ip'ip production from such configurations. 



A. Inycccc) Configurations 

In this section, we follow the arguments of Ref. The J/tp production cross section from an \nvcc) configuration 

is 

atihN) = f^/f,ai,{hN) , (25) 

where is the fraction of the intrinsic cc pairs with mass below the DD threshold that coalesce into a final-state 
J/ip, cc — > J/tf}. Recently we estimated /^/^ ~ 0.03 and f^fp w 0.014 ^ with nic — 1.5 GeV. The predicted cross 

sections are af^{Tr~N) = 15 nb and af^ijpN) = 9.8 nb, in agreement with the NAB single J/ip data Given the 
intrinsic charm cross section for single J /ip production, we use the ipij) data to estimate the probability to produce 
two cc pairs in the projectile, Pjcc, assuming that all the ijjip pair events arise from this mechanism, as in Ref. 
Then 



'(^^) = /|/.^^ic(/*^) = U/H^<{hN) . (26) 

J 1 r. J i r 



The measured value, cr^^ = afJ'{Tr^N) w 20 pb [|l|, requires Pice ~ 4.4%Pic, implying that it is relatively easy 
to produce a second intrinsic QQ state if one is already present. If Pec is independent of the projectile identity, 
afJ'{pN) Ki 6.1 pb, only 20% of the measured cross section at 400 GeV Note that if part of the 7t~N tpip 
production cross section is assumed to arise from leading-twist mechanisms. Pec is reduced. 
The single J distribution in this configuration is 

— — / dxidx ~] '^(^^i Xq-^ Xc^)(5(x^,2 Xc2 Xc2^ • (^'^) 

^^^01 J T uXi . . . UXji 
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We find {x^) = 0.36 for the pion projectile and {x^) = 0.33 for the proton. The number of single J/tp's is twice 
the number of ipip pairs. The single J/t/i's in the 'lp^p pairs have a lower average momentum fraction than those from 
\nvcc) Fock states, where (x.^) — 0.62 for a pion and (x^) — 0.51 for a proton [^,^. The V'V' Pa^ir distribution is 

dx^^- o{x^4,^ - x^^ - x^^) . (28) 



dXqpip J dxip^ 

Here we find (x^^) — 0.72 for the pion and (x^^) — 0.64 for the proton. We compare the frame-independent 
calculation to the combined 150 and 280 GeV data. The ijji! pair distributions are in Figs. 10(a) and 11(a) while the 
single J/tp distributions are shown in Figs. 10(b) and 11(b). 
The Tpif) pair mass distribution from this configuration is 

= / n '^^^'^^^T., n ^dmlJ^kT,^,, ^^d^kr,^^ "S'i/ZA-^ ^^^^ 

(m^^. m'^c '^Tc\ ^ ^ ^ 

X S — — S{kT.cj + kr^c, - kT.,pj)5{x^^ - Xc^ - Xc^) 

\ ^ipj ^Cj ^Cj J 

c ( ^T,ipip "^T.Vi '^T,^/.2 1 r/r , r r ^^^^ \ 

X 0(kT,^i + Kt,-!/J2 - kT,4,4,)0{x^.^ - X^^ ~ X^2) I 

\ X^ip Xip-^ Xil;2 I 

where 2mc < < 2m]j. The delta functions insure conservation of momentum for the single J/^ constituents and 
the pair. The results are compared with the data in Figs. 10(c) and 11(c). We find {M.^^) m 7.7 GeV for the pion 
beam and {M^i,^) « 7.4 GeV for the proton beam. The smaller {M^,^) in pN interactions is compatible with the data 
01- 

This mechanism provides a simple frame-independent source of fast ipip pairs as well as fast single J/V''s and 
reproduces the strongly correlated behavior of the pairs. The general agreement with the intrinsic charm model 
is good. 



B. \nvccbb) Configurations 

We now discuss ipip production from \nvccbb) configurations. A projectile in such a configuration could produce a 
ip'ip PE^ir from a combination of i? ^ J/t/j decays with cc — *■ J/ip coalescence. The ipip production cross section is 

<^ttihN) = U/hj^aiUhN)BR{B ^ J/^jX) . (30) 

If we assume Pkb/^'ib ~ Picc/Pic, then o'j'^b (tt^ A^) « 0.044 pb and crfctiP^) ~ 0-020 pb. These cross sections are 
equivalent to those from leading-twist bb production, much less than the measured cross sections. The wN cross 
section is larger since f^/p < /^/tt. 

The J/V' distribution from this state is 

^ = / n d^^d^^j^"^^ - - (31) 

^ i—l 

since the intrinsic cc coalescence is unaffected by the b quark hadronization. We find {x^) — 0.31 from a pion projectile 
and (a;^) = 0.28 from a proton, independent of frame. These J/V' distributions are shown in the dot-dashed curves 
in Fig. 12(a) and 13(a). 

If the B meson is produced by fragmentation, 



dp 



icb 



dx 



Bf 



i=l 

Alternatively, if the b quark hadronizes by coalescence. 



/^^ dB' 
dzY\ dxi- — — dx^S{x^ - Xc - x-)DB/biz)S{xBp - zxb) ■ (32) 
dxi...dxn 



Y[dxi- ^- — dx^S{x^ - Xc- Xc)S{xbc - Xb- xi) . (33) 



dPcb /■A , dPi 

dxBc 
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Again, note that the B production is frame independent but the J/tp distributions are presented in the laboratory 
frame. The Bp and Be distributions are given in Figs. 12(a) and 13(a) for the pion and proton while the J/^p 
distributions from the decays appear in Figs. 12(b) and 13(b). 

The ipip pairs from the \nvccbb) Fock state arise from a B ^ J/ipX decay combined with cc — > J/ip coalescence 
and boosted to the laboratory frame. The average ipip momentum is larger when the parent B was produced by 
coalescence. The results are shown in Figs. 12(c) and 13(c). The average momentum fractions of the B mesons and 
the iptp pairs from this configuration are given in Table III. While the average longitudinal momentum fraction is 
relatively large, the pair mass is sharply peaked near threshold with {M^p^) k, 6.5 GeV, as shown in Figs. 12(d) and 
13(d). Additionally, the cross section is a factor of 500 below the data. Therefore it is unlikely that this mechanism 
could produce the -^pij) pairs observed by NA3. 



3. Leading-Twist J /ip Production with an Intrinsic Charm Configuration 

We have also studied the possibility of double J/ ij] production by simultaneous leading-twist fusion and coalescence 
of an intrinsic cc pair, cc — > J /ip. In this case, the parton distributions of the projectile pion and proton are obtained 
from the \ncc) Fock configurations, as shown in Fig. 1(g), where n represents the valence quarks plus other light 
partons. Thus for qq annihilation with a pion valence quark, only the lowest intrinsic charm Fock state, judcc), is 
needed. For gluon fusion, an additional gluon is necessary, \udgcc). Annihilation with a sea quark requires at least 
a six particle Fock state, \udqqcc). With a proton projectile, the lowest-lying configurations have five, six, and seven 
partons. The valence quark probability distribution from an intrinsic cc — > J /ip state is given by 

qv„{xi) = '^-^ = [ dx2pdxcdx-dx^2- — — 6{x^^ - Xc- x-) , (34) 

axi J dxi . . . dxn 

where X2p is a light quark in the projectile (the 2p subscript is used to distinguish it from the target momentum 
fraction) and Xi is the projectile parton contributing to the fusion process, either a valence or sea quark or a gluon. 
We also need the "02 distribution when xi is fixed by the fusion process, 

qv^ixiix^^) = —^^^ — 1^1= / dx2pdxcdxc , '^^^"j — 5{x^^ - Xc- Xc) . (35) 
axiax02 J axi . . . axn 

The sea quark and gluon distributions are similarly calculated. The target parton distribution functions are GRV LO 
pof and Duke-Owens 1.1 | [T^ . The projectile and target momentum fractions in the leading- twist fusion process are 

xi.2 ~ (±a;^i + ^J'x^~{~4t^)/2, where Imd \fs < r < 2mr)/s/s. The cc pair is assumed to neutralize its color and 

produce a J /i}) by emitting a soft gluon without significant momentum loss. 
The J /tp produced by parton fusion has the momentum distribution 

da /•2mc/yj 

-1 = '^Ff^/h / -^=^=[qvp{x2)qv^{xi)aqq^cc + Gp{x2)GT,{xi)agg^cc 

U'XtPx J2me/\/i X2Jx'^^^+At'^ 

+qs^{x2)qs^{xi)aqq^cc] , (36) 

written in order of increasing number of partons in the minimal projectile Fock state. We define xq^^i^x) — xup(x) + 
xdp{x) for annihilation with the ■k~ iud) valence quarks and xqsp{x) = xup{x) + xup(x)+xdp{x) + xdp{x) + 2xsp{x) for 
annihilation with pion sea quarks where xup = +xus and x dp = xd^ +xds . These J/V' distributions, shown 
in Figs. 14(a) and 15(a) for pion and proton beams, are central. In the laboratory frame, {x^)p{ « 0.26 from n^p 
interactions and {x^)p{ « 0.18 from pp interactions, nearly independent of the beam energy. The J/tp distribution 
from the intrinsic charm state is 

da „ i'2mo/V^ 

-■[qvpix2)qv„{xi;x^^)agg^cc 



= '2Ff^/h / dx^^ / 



I Jtp/!L I (/^l I / 

Jo J2mJ^ X2Jx^^+'±, 

+ Gp{x2)G 

+ qsp {X2) ] ■ (37) 



These distributions, shown in the laboratory frame in Figs. 14(b) and 15(b), have larger average momentum fractions, 
(a;^)ic ~ 0.56 and 0.46 for pion and proton projectiles. Thus one fast and one slow J/ip are produced, in agreement 
with the single J/ip data. The tptp distributions are strongly forward peaked, with (x^^) « 0.84 for a pion and, for a 
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proton, {Xjp^) « 0.68, larger than the data. The average ipi' P^^ir niass is « 8 GeV. However, the mass distributions 
do not appear strongly correlated. 

The parameter F is introduced since the leading-twist fusion cross section includes all cc resonances below the DD 
threshold. When rric — 1-5, F « 1/5. We also include the factor /^//j for J/ip production from the intrinsic charm 
state. These two suppression factors predict a iptj) production cross section of « 3.5 pb for 7rp and « 2.5 pb for pp 
production. 



IV. Summary 

A summary of all the predicted cross sections is given in Table IV. From the magnitude of the calculated cross 
sections alone, it appears that approximately 50% of the measured ^tl) cross section can be attributed to the leading- 
twist mechanisms. However, only the higher-twist mechanisms produce fast ipij: pairs. The intrinsic heavy quark states 
offer a promising alternative, especially since all pairs produced through these mechanisms carry a significant 
fraction of the projectile momentum. Particularly, the \nvcccc) configuration can also account for the size of the ipij} 
cross section, as discussed in ||]. 

The correlations in BB production also suggest an intriguing test of the intrinsic heavy quark mechanism. In 
particular, if both B's are produced by coalescence, all the pion momentum would be transferred to the BB pair. The 
same type of correlation should be observable in DD production in nA interactions. Measurements of the D~ /D~^ 
asymmetry could be extended to the momentum and invariant mass distributions of DD pairs. 

We thank S.J. Brodsky, P. Hoyer, G. Ingelman, and W.-K. Tang for stimulating discussions and J.-C. Peng for 
encouragement and for providing the B decay Monte Carlo. 
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Figure Captions 



Figure 1. (a) Examples of 0{a'^) QCD graphs for gg tjjijj. (b) The diagrams for 0{a'^) qq tjjijj production, 
(c) Leading-twist bb production by qq annihilation and gg fusion, (d) Drell-Yan type resonance production. The 
intermediate state contains four quarks and ip denotes a color-octet-vector cc. (e) An intrinsic bb pair illustrating 
the coalescence {Be) and fragmentation {Bp) production processes, (f) An example of a pair of intrinsic QQ states 
in the projectile, (g) Diagrams of leading- twist J/ip production with the projectile in an intrinsic cc state. Fusion 
production by valence quark annihilation, gluon fusion, and sea quark annihilation are shown respectively. 

Figure 2. The 0{a'^) 'K~p iptp distributions in the laboratory frame compared to the NA3 data. The x^p^p 
distributions are shown in (a), the distributions in (b), and the pair mass distribution in (c). The curves show: at 
150 GeV, GRV LO (solid) and DO 1.1 (dashed); at 280 GeV, GRV LO (dot-dashed) and DO 1.1 (dotted). 

Figure 3. The 0{a'^) pp tpip distributions in the laboratory frame. The x^^ distributions are shown in (a), the 
distributions in (b), and the pair mass distribution is compared to the NA3 data in (c). The curves show: at 400 
GeV, GRV LO (solid) and DO 1.1 (dashed); at 800 GeV, GRV LO (dot-dashed) and DO 1.1 (dotted). 

Figure 4. Loading-twist ir^p bb prochiction in the laboratory frame, compared to the NA3 ^jiJj data. The .r^^^ 
distributions are shown in (a) and the x^ distributions in (c) with GRV LO (solid) and DO 1.1 (dashed) at 150 GeV 
and GRV LO (dot-dashed) and DO 1.1 (dotted) at 280 GeV. The remaining results, with the GRV set, compares the 
assumptions of bb decay vs. BB production and decay. The solid and dot-dashed curves give the bb distributions at 
150 and 280 GeV while the dashed and dotted curves show the BB distributions at the same energies. The x^^ and 
x^ distributions are shown in (b) and (d), the tjjil) pair mass distribution in (f), and the invariant mass of the bb and 
BB pairs in (e). 

Figure 5. Leading- twist pp bb production in laboratory frame. The x^-^ distributions are shown in (a), the Xb 
distributions in (c) with GRV LO (solid) and DO 1.1 (dashed) at 400 GeV and GRV LO (dot-dashed) and DO 
1.1 (dotted) at 800 GeV. The remaining results, with the GRV set, compares the assumptions of bb decay vs. BB 
production and decay. The solid and dot-dashed curves give the bb distributions at 400 and 800 GeV while the dashed 
and dotted curves show the BB distributions at the same energies. The x^^ and x^ distributions are shown in (b) 
and (d), the iptp pair mass distribution in (f), compared to the NA3 data, and the invariant mass of the bb and BB 
pairs in (e). 

Figure 6. The 'jt~N — > 2~^~^cccc tp^p distributions in the laboratory frame, compared to the NA3 data. The x^^ 
distributions are shown in (a), the x^ distributions in (b), and the M^^ distribution in (c). The calculations are: at 
150 GoV, GRV LO (solid) and scaling (dashed) gluon distributions; at 280 GeV, GRV LO (dot-dashed) and scaling 
(dotted) gluons. 

Figure 7. The pN 2'^~^cccc ipip distributions in the laboratory frame. The x^^ distributions are shown in (a), 
the x^ distributions in (b), and the M^^ distribution is compared with the NA3 data in (c). The calculations are: at 
400 GeV, GRV LO (solid) and scaling (dashed) gluon distributions; at 800 GeV, GRV LO (dot-dashed) and scaling 
(dotted) gluons. 

Figure 8. Frame-independent bb production is calculated from an intrinsic \udbb) configuration. In (a), the x^^-g^ 
(solid) and x^^^^ (dashed) distributions are shown. If both S's are produced by coalescence, x^^-g^ = 1. In (c), 
the XBp (solid) and xbc (dashed) distributions are shown. The invariant masses of the bb (solid), BpBp (dashed), 
Be Bp (dot-dashed), and Be Be (dotted) pairs are shown in (e). The J ftp distributions are in the laboratory frame. 
The distributions are shown in (b). The solid and dot-dashed curves show BpBe tpip pairs at 150 and 280 
GeV while the dashed and dotted curves show BpBp ijjtjj pairs at 150 and 280 GeV. The thin solid and dashed 
curves are results from BeBe V'V' decays at 150 and 280 GeV. The x^ distributions are in (d). The Be J/tpX 
decays are given in the dashed (150 GeV) and dot-dashed (280 GeV) while Bp — > J/ipX decays are shown in the 
solid (150 GeV) and dotted (280 GeV) curves. Finally, the M^^ distributions from BpBp decays (solid) and BcBc 
decays (dashed) are given in (f). 
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Figure 9. Frame-independent bb production is calculated from an intrinsic \uudbb) configuration. In (a), the x^^-g^ 
(solid) and x^^^^ (dashed) distributions are shown. In (c), the XBp (solid) and xbc (dashed) distributions are shown. 

The invariant masses of the bb (solid), BpBp (dashed), and BcBp (dot-dashed) pairs are shown in (e). The J/ij) 
results are shown in the laboratory frame. The x^^p distributions are shown in (b). The solid and dot-dashed curves 
show the BpBc V'V' pairs at 400 and 800 GeV while the dashed and dotted curves show BpBp —* tptp decays at 
400 and 800 GeV. The distributions are in (d). The Be — > J/tpX decays arc given in the dashed (400 GeV) and 
dot-dashed (800 GcV) and the ^ J/ipX decays are given by the soUd (400 GeV) and dotted (800 GeV) curves. 
Finally, the M^^ distribution is in (f). 

Figure 10. Frame- independent tptp production from a \udcccc) configuration. The x^^ distribution is shown in (a), 
the x^ distribution in (b), and the mass distribution in (c). 

Figure 11. Frame- independent ipijj production from a \uudcccc) configuration. The x^^ distribution is shown in (a), 
the x^ distribution in (b), and the mass distribution in (c). 

Figure 12. Results for J/ip and tpip production from an intrinsic \udccbb) configuration. In (a) we show the frame- 
independent results for: cc J/ij.) (solid). Be (dashed) and Bp (dot dashed) production. All other results arc in 
the laboratory frame. In (b), the Be decays are given by the solid and dotted curves for 150 and 280 GeV while the 
Bp decays are shown in the dashed and dot-dashed curves for 150 and 280 GeV. In (c) and (d) the a;^^ and M^^ 
distributions result from the combination of cc J/t/i with a decaying B: pairing with a Be results in the dashed 
(150 GeV) and dotted (280 GeV) curves while a Bp pairing produces the results given in the solid (150 GeV) and 
dot-dashed (280 GeV) curves. 

Figure 13. Results for J/^/j and production from an intrinsic \uudccbb) configuration. In (a) we show the frame- 
independent results for: cc J/ip (solid), Be (dashed) and Bp (dot dashed) production. All other results are in 
the laboratory frame. In (b), the Be decays are given by the solid and dotted curves for 150 and 280 GcV while the 
Bp decays are shown in the dashed and dot-dashed curves for 150 and 280 GeV. In (c) and (d) the x-^.^ and M-^^ 
distributions result from the combination of cc J /ip with a decaying B: pairing with a Be results in the dashed 
(150 GeV) and dotted (280 GeV) curves while a Bp pairing produces the results given in the solid (150 GeV) and 
dot-dashed (280 GeV) curves. 

Figure 14. ipxp pairs produced by simultaneous J/V' production by leading-twist fusion with a pion in an judncc) 
configuration, in the laboratory frame. The leading- twist and intrinsic cc J/V' distributions are given in (a) and (b), 
and the x.^^ and mass distributions are shown in (c) and (d). The calculations are: at 150 GeV, GRV LO (solid) and 
DO 1.1 (dashed); at 280 GeV, GRV LO (dot-dashed) and DO 1.1 (dotted). 

Figure 15. -i/'V' pairs produced by simultaneous J /ip production by leading- twist fusion with a proton in an \uudncc) 
configuration, in the laboratory frame. The leading- twist and intrinsic cc J/ip distributions are given in (a) and (b), 
and the a;^^ and mass distributions are shown in (c) and (d). The calculations are: at 400 GeV, GRV LO (solid) and 
DO 1.1 (dashed); at 800 GeV, GRV LO (dot-dashed) and DO 1.1 (dotted). 
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TABLE I. The average laboratory momentum fractions and invariant masses for leading-twist 
BB production and decay in hp interactions calculated with the GRV LO parton distributions. 

Note that the combined 7r~A'" ipipX data give (a;^) = 0.33, (x^^) = 0.66, and (M^^) = 7.4 GeV 
while the average tpi^ mass from the pN tptpX measurements is (M^^) = 6.8 GeV. 



h, Pbeam (GeV) 


(xb) 








(M^^) (GeV) 


(M^V.) (GeV) 


TT", 150 


0.34 


0.24 


0.65 


0.50 


11.5 


7.1 


TT", 280 


0.27 


0.19 


0.52 


0.41 


12.1 


7.4 


p, 400 


0.19 


0.13 


0.37 


0.29 


11.9 


7.3 


p, 800 


0.15 


0.10 


0.27 


0.22 


12.3 


7.5 



TABLE II. The average momentum fractions of B and BB production and decay from intrinsic \nvbb) states in the projectile 

h. The B results arc frame independent, the and 1/;^" averages arc in the laboratory frame. The invariant masses arc given 
in GeV. Recall that Bf denotes B mesons produced by fragmentation while Be denotes mesons produced by coalescence with 
projectile valence quarks. Note that the combined 7r~iV —f iptpX data give {x^) = 0.33, (xy,^) = 0.66, and (My,^) = 7.4 GeV 
while the average iptp mass from the pN — > iptpX measurements is (M^^) = 6.8 GeV. 



Bp Be BpBp BpBc BcBc 



h, Pbeam (GeV) 


(xb) 




{xb) 


{x^} 


{■'■bb) 


{^^Iplp ) 


{^bb) 


(^^tplp } 


{^bb) 




TV', 150 


0.33 


0.35 


0.50 


0.44 


0.64 


0.71 


0.80 


0.76 


1.0 


0.81 


TT", 280 


0.33 


0.31 


0.50 


0.40 


0.64 


0.60 


0.80 


0.67 


1.0 


0.75 


p, 400 


0.29 


0.26 


0.43 


0.35 


0.55 


0.51 


0.69 


0.58 






p, 800 


0.29 


0.24 


0.43 


0.32 


0.55 


0.44 


0.69 


0.52 
















(Mbb) 












7r~ 










12.5 


7.65 


12.2 


7.5 


11.8 


7.05 


P 










12.5 


7.65 


12.4 


7.6 







TABLE in. The average momentum fractions for B production and decay from an intrinsic cchb state. The ipip pairs 
are produced through cc — > J/xp coalescence and B — > J/xpX decay. The B averages are frame independent, the J/il) and 
tpil) averages are in the laboratory frame. Recall that Bp denotes B mesons produced by fragmentation while Be denotes 
mesons produced by coalescence with projectile valence quarks. Note that the combined 7r~7V ^ iIj'^X data give (xy,) = 0.33, 
{x^^) = 0.66. 



h, Pbeam (GeV) 


{xb) 


Bf 


{x^) 


{xb) 


Be 


{Xi,) 


ipBp 

{Xipii) 




TT", 150 


0.23 




0.30 


0.35 




0.36 


0.69 


0.74 


TT", 280 


0.23 




0.25 


0.35 




0.32 


0.62 


0.67 


p, 400 


0.21 




0.22 


0.32 




0.28 


0.54 


0.60 


p, 800 


0.21 




0.19 


0.32 




0.25 


0.49 


0.56 



TABLE IV. A summary of a^,p for all the studied processes. When only one value is given for 
more than one energy, the results arc nearly energy independent. 



Process 


0,1, 


,,(^-7Y) (pb) 


(y^iAv^) (pb) 






150 GeV 


280 GeV 


400 GeV 


800 GeV 


NA3 


18 ±8 


30± 10 


27± 10 




0{at) 
0{al) bb 


7.4 


16.5 


7.2 


26.0 


0.028 


0.24 


0.03 


0.16 


2^^cccc 


1.7 


10.0 


12.8 


55.2 


\nvbb) 




0.36 




0.51 


\nvcccc) 




s20 




6.1 


\nvccbb) 




0.044 




0.020 


ncr) 






1.9 


2.8 
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